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Abstract 

A measurement of correlations between event-plane angles <t>„ is presented as a function of centrality for Pb- 
Pb collisions at = 2.76 TeV. These correlations are estimated from observed event-plane angles *¥„ obtained 
from charged particle or transverse energy flow measured over a large pseudorapidity range \tj\ < 4.8, followed by a 
resolution correction that accounts for the dispersion of ¥„ relative to <$>„. Various correlators involving two or three 
event planes with acceptable resolution are measured. Significant positive correlations are observed for 4{<&2 - 
6(<E>2 - 4*6), 6(<1>3 - 6 ), 2®2 + 3O3 - 50 5 , 20 2 + 4®4 - 6<J> 6 and -10O 2 + 40 4 + 60 6 . However, the measured 
correlations for 24>2 - 6O3 + 4C>4 are negative. These results may shed light on the patterns of the fluctuation of the 
created matter in the initial state as well as the subsequent hydrodynamic evolution. 
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In recent years, the measurement of harmonic flow coefficients v n has provided important insight into the hot 
and dense matter created in heavy ion collisions at the Relativistic Heavy Ion Collider (RHIC) and the Large Hadron 
Collider (LHC). These coefficients are generally obtained from a Fourier expansion of particle production in azimuthal 
angle <p, ^ oc 1 + 2J]" =1 v„ cos n(cf> - <£>„), where <E>„, known as the event plane (EP), represents the phase of v„. 
Experimentally, O,, is estimated using the azimuthal distribution of particles in the event. The estimated EP direction 
T*,, is referred to as the observed event plane, to distinguish it from <£„. 

Besides the v 2 coefficient, which reflects mainly the hydrodynamic response of the matter to the elliptic shape of 
the overlap region, significant first order (vi) and higher order (V3-V6) coefficients were also obtained EIUE]. These 
non-elliptic harmonics suggest that the matter densities can fluctuate in the initial state, and the harmonics simply 
reflects the hydrodynamic responses to these fluctuations. More details about these fluctuations can be inferred from 
the correlations between 0„ of different order. Calculations based on Glauber model suggest the existence of signif- 
icant correlations in the initial state [4, 5, 6, 7]; the correlations can also be generated dynamically in hydrodynamic 
evolution due to non-linear effects [8 1. This proceedings present measurement of EP correlations |9|, which can shed 
light on these initial state fluctuations and dynamic mixing between different <£„ in the final state. 

The correlation between several EPs are characterised by a set of cosine functions, (cos(ci<E>i + ... + lc/(S)i)) with 
c\ + 2c2--- + lc[ = O. Two-plane correlators can be written as (cos A:(<t>„ - <£,„)), where k is the least common multiple 
of n and m. Experimentally, these correlations are obtained from the observed event planes as [7|: 

(cosfe-S, + ... + 1cm) = t°1 C 'I\ + p V Q ;m^ Res{c„«>P„} = (cos c^P, - <D„)> (1) 

Res{ciTi)...Res{c / /*I'/} 
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where the resolution factors Res{c„77 v F„} are determined using the standard two-subevent (2SE) or three-subevent 
(3SE) methods [10]. To avoid auto-correlations, the needs to be measured using subevents covering different 77 
ranges, preferably with a gap in between. A large number of correlators could be studied. However, the measurability 
of these correlators is dictated by the value of ResfjnT,,}; they are measurable in ATLAS for n = 2-6 and j values 
up to 6 [9]. This analysis focuses on two- and three-plane correlations. The latter can be generally constructed from 
several two-plane relative angles, and hence probes the correlations of two angles relative to the third Q, such as: 



2<I> 2 + 44>4 - 6®6 = 4(<&4 - 2 ) - 6(<&6 - $2), — 10<t> 2 + 40 4 + 6<J> 6 = 4(® 4 - <E> 2 ) + 6 (3>6 
2<sin4(0 4 - 2 ) sin6(<l>6 - 2 )> = <cos(2<t> 2 + 4$ 4 - 60 6 )> - (cos(-10a> 2 + 4<t> 4 + 60 6 )> 



<*> 2 ) 



(2) 



In a two-plane correlation (2PC), the event is divided into two subevents symmetrically around r\ — with a gap 
in between such that nominally they have the same resolution. Each subevent provides its own measurement of the 
two event planes: and for positive r\ and W% and ^ for negative 77. This leads to two statistically independent 
measurements of the correlator, which are averaged to obtain the final signal: 



<cosfc(O n -<I> m )> = 



(cos k(V p n - + (cos k(V N n - *0) 
Res{^P}Res{A:^} + Res{/t*I'N}Res{/fcT' p , 



(3) 



In a three-plane correlation (3PC), three non-overlapping subevents, labeled as A, B and C, are chosen to have 
similar 77 coverage. In this analysis, subevent A and C are chosen to be symmetric about 77 = 0, and hence have 
identical resolution, while the resolution of B can be different. The symmetry between A and C reduces the 6 inde- 
pendent combinations into three pairs of measurements, labeled as Typel, Type2 and Type3. For example, the Typel 
measurement of the correlation 2<5 2 + 3<t>3 - 5*1)5 is obtained by requiring the V F 2 plane to be measured by subevent B: 



<cos(20 2 + 3<D 3 - 5<6 5 )> Typel = 



(cos(2*Pf + - 5^)) + (cos^f + 3Tf - 5Wj)) 
Res{2 ( Tf}Res{3 > I'^}Res{5 ( I'^} + Res{2 > Tf }Res{3^}Res{5^} 



(4) 



Similarly, the Type2 (Type3) measurement is obtained by requiring the ^3 (^s) to be measured by subevent B. 

This analysis is based on approximately 8 /jb -1 of minimum bias Pb-Pb data collected in 2010 at ~^s^ - 
2.76 TeV 0. Three subsystems of ATLAS IfTTII are used to measure the event plane: the inner detector (ID) over 
I77] < 2.5, the barrel and endcap electromagnetic calorimeters (ECal) covering I77I < 3.2 and the forward calorimeter 
(FCal), which extend the calorimeter coverage to I77I < 4.9. These detectors are divided into a set of segments each 
covering 0.5 to 0.7 units in 77, and the subevents are constructed by combining these segments. The subevent combi- 
nations for the two- and three-plane correlations are listed in Table [Tj The combination of ECal and FCal (ECalFCal) 
provides the default measurement, while the ID serves as cross-checks. 



Two-plane correlation (A,B) 


Default 


ECalFCal^ 77 €(0.5,4.8) 


ECalFCaliv 77 £(-4.8,-0.5) 


Cross-check 


JD P 77 e(0.5,2.5) 


ID N 77 £(-2.5,-0.5) 


Three-plane correlation (A,B,C) 


Default 


ECalp 77 6(0.5,2.7) 


FCal 6(3.3,4.8) 


ECaL, 77 £(-2.7,-0.5) 


Cross-check 


IDp 77 6(1.5,2.5) 


ID 77 e(- 1.0, 1.0) 


IDjv 77 £(-2.5,-1.5) 



Table 1 : The definition of the subevents and their i] coverage using the two-plane and three-plane correlations. 

Figure [T] shows the two-plane relative angle distributions for the (20-30)% centrality interval from the FCalECal. 
To estimate the detector effects, a background distribution is calculated by combining (or W,„) in one event with 
T*,,, (or *P„) in another event with similar centrality (matched within 5%) and primary vertex z-position (matched 
within 3 cm). The correlation function is then obtained as the ratio of same event correlation (S) to the background 
distribution (B): 



C(kQ¥ n - = 



S(fcQF„ -Y m )) 
B{kQ¥ n - T m )) 



The correlation functions show significant positive signals for 4( V F 2 - *I / 4 ), 6( V F 2 - ¥3), 6( V F 2 - ^d) and 6(^3 
The observed correlation signals are calculated as the cosine averages of these correlation functions. 



(5) 
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Figure 1: Relative angle distributions between two event planes for the (20-30)% centrality interval for foreground (open circles), background 
(open squares) and correlation function (filled circles). The correlation functions are used via Eq.[3]to obtain the two-plane correlators. 



The resolution factors Resjy'nT,,} are determined using the 2SE and 3SE methods, with their differences as part of 
the systematic uncertainty. The potential influence of non-flow effects due to short-range correlations are studied by 
requiring a minimal 77 gap, 77 m ; n , between the two ECalFCal subevents. Results obtained for 77 m i n values in the range 
of to 6 are found to be consistent. In all cases, the influences of these short-range correlations are negligible for 
?7min > 0.4 (the choices of the subevents in Table[T]have a minimum 77 gap of 0.6). 

The analysis procedure discussed above is also valid for the 3PC analysis, except that the three independent mea- 
surements of 3PC for each correlator need to be combined. Figure|2]shows the relative angle distributions for various 
three-plane correlators from the Typel measurement in the (20-30)% centrality interval. The observed correlation 
signals are calculated as cosine average of the correlation function in a generalization of Eq. [5] 



C(c„n < F„ + c m m^ m + e,,/!^) = 



SjCnll^n + C m m x ¥ m + C;,/W;,) 

B(c n riY n + c m mW m + c h h^ h ) 



(6) 



where the background distribution is constructed by requiring and from different events. The correlation 

functions show significant positive signals for 2 X F 2 + 3^ - 5^5, 2 V P 2 + 4*P 4 - &¥ 6 , and -10*F 2 + 4 V P 4 + 6 X ¥ 6 , while 
the signal for 2 v t , 2 - + 4^4 is negative, and the signals for the remaining correlators are consistent with zero. 

Figures [3] and |4] show the centrality dependence, given here by average number of participating nucleons (N part ), 
of the two-plane and three-plane correlators, respectively. The systematic uncertainties are strongly correlated across 
the centrality range shown. Given that the ID is a different type of detector and measures only charged particles, 
the consistency between the ID and ECalFCal attests to the robustness of this measurement. Large positive values 
are observed in most cases and their magnitudes usually decrease with (-/V part ). The value of (cos 6(<I>2 _ $3)) is 
small (< 0.02), yet exhibits a similar dependence on (N pm -i). Interestingly, two other correlators show very different 
trends: the value of (cos 6(03 - fl>6)) increases with (N pai t), and the value of (cos(2<l>2 - 6O3 + 4<J>4)) is negative 
and its magnitudes decreases with (A^ par t)- The centrality dependence trends for some correlators are qualitatively 
similar to the predictions of the Glauber model 0161121, such as (cos 4(0 2 - O4)}, (cos 6(0 2 - 6 )>, < cos 6(<I>2 - <& 3 )>, 
(cos(2<I>2 + 3<I>3 - 5<t> 5 )) and (cos(2<L>2 + 4<I>4 - 60s)); however, the Glauber model fails to reproduce the magnitudes 
and centrality dependences for other correlators such as (cos 6(<I>3 - <t>6)), (cos(2i>2 - 6O3 + 4O4)) and 
(cos(-10<l>2 + 4<I>4 + 6<t>6)). This behavior might suggest the importance of non-linear response of the hydrodynamic 
evolution of the medium to the fluctuations in the initial geometry (8). 

In summary, measurements of several correlators between two and three event planes are presented using 8 /ib~' 
Pb-Pb collision data at *^s^ = 2.76 TeV f9]. Several different trends in the centrality dependence of these correlators 
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Figure 2: Relative angle distributions between three event planes for the (20-30)% centrality interval for Typel correlation for foreground (open 
circles), background (open squares) and correlation function (filled circles). The correlation functions are used via equations similar to Eq. [4]to 
obtain the three-plane correlators. 
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Figure 3: The final corrected eight two-plane correlators (cos jk(<t> n - <D„,)) as a function of (iVpart) measured by ECalFCal and ID. The middle 
two panels in the top row have j = 2 and j = 3, respectively, while all other panels have 7 =1. The error bars and shaded bands indicate the 
statistical uncertainty and total systematic uncertainty, respectively. 



are observed. Some of the trends are qualitatively similar with predictions based on the Glauber model, while others 
differ significantly. These observations suggest that both the fluctuations in the initial geometry and dynamical evo- 
lution of the medium in the final state are important for creating these correlations in momentum space. A detailed 
theoretical description of these correlations may lead to new insights into the patterns of the fluctuation of the created 
matter in the initial state as well as the subsequent hydrodynamic evolution in heavy ion collisions. 
This work is in part supported by NSF under award number PHY-1019387. 
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Figure 4: The final corrected six three-plane correlators as a function of (Np m t) measured by ECalFCal and ID. The error bars and shaded bands 
indicate the statistical uncertainty and total systematic uncertainty, respectively. 
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